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TITLE 

TYPE 

ACCURACY 

DURATION 

METHOD OF USE 



Automatic Inversion of a Symmetric Matrix 
Entire program 

Depends on the condition of the matrix to be inverted 
See text 

The code is read into the memory in the usual way. After 
a short time, Illiac will stop. Read in the elements 
of the matrix and parameters . Replace original tape 
in reader and restart computer. Since A is symmetric, 
its inverse is symmetric, and it will "be necessary to 
print only half the off-diagonal elements. The first 
number that appears is the row number. This is followed 
by the diagonal elements and the upper off -diagonal 
elements. The inverse of the matrix should be scaled 
by the same factor by which the original has been scaled. 
A sample of the final result follows: 



PUNCHING OF DATA 



7^95*858^9 
lla75«086"96 
18^55.762^7 

15562.1t6767 

Since the original matrix was scaled down by ten, 
the above elements must be scaled down by ten. 
The elements should be scaled such that the sum of 
the squares of the matrix elements is less than one 
fourth. Since A is a symmetric matrix, it is necessary 
to store only half the off-diagonal elements . The 
lower off -diagonal elements and diagonal elements are 
punched, row by row, as a sign followed by up to 12 
decimal digits. The final fraction in the sequence 
is terminated by a J. This is followed by fF where 
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f determines the number of places of the fraction 
desired to be printed. Finally, nN appears on the 
tape where n is the order of the matrix. The data 
tape should appear as follows: 

a 00 

a 10 * a ll 

a n _ , *.., a t . JfFnN 
n-1, n-1, n-1 

DESCRIPTION The symmetric matrix A is diagonalized "by a series 

of orthogonal transformations . The final matrix 



product can "be 


written 


D = 


B T AB. 


Then 

A = 


m „1 rn 

B DB (BB 


= 


T 
BDB 


and 




A' 1 


„-T -1 -1 
= B D B 




-1 T 
= B D B . 



= I) 



Since D is a diagonal matrix, the inverse has for 
elements the reciprocals of the diagonal elements of D. 
After Routine k2 finds the matrices B and D, this 
routine inverts D, performs the necessary matrix multi-* 
plication, and then prints out the result. This routine 
is essentially a combination of routines X 1, 21, k2, 
P 2, P k, N 2» It is also possible to invert an asym- 
metric matrix using this routine. A matrix C can he 
made symmetric by multiplying it by its transpose. 
Thus 



Then 



and 



T 
CC = A s where A is symmetric 



a" 1 = cf T c" 1 CD 



C" 1 = C T A" 1 . (2) 



Operations (l ) and (2 ) can be performed by use of 
Routine Ml - 91. 



-3- 



LENGTH OF TIME 



CONDITION OF THE MATRIX 



The length of time of the program, the accuracy, 
and the maximum matrix allowed to he inverted are 
essentially dependent on Routine 42. The maximum 

matrix that can "be handled is of order 23 • Routine 

3 -3 
k-2 takes approximately 5n x 10 seconds to perform 

one iteration. The number of iterations varies directly 

•with the order of the matrix , from about four iterations 

for a third order matrix to 7 for a twenty-third 

order matrix. To perform its various operations and 

print out the results, this routine requires approximately 

kn (n + l) (3n + 63 + 5f ) milliseconds where n is the 

order of the matrix and f is the number of places the 

fraction is printed- It was found that slightly more 

than Ik minutes were necessary to invsrt a matrix of 

order 23 and print out the result. 

The accuracy of the result will naturally depend on 

the conditioning. The conditioning of the matrix is 

defined here as the ratio "between the smallest and 

largest eigenvalues. As this ratio becomes very small, 

the accuracy becomes very poor. This ratio is printed 

as the last number and is followed by an N« The number 

of figures of accuracy will be no greater than 11 

minus the number of zeros that appear in the rounded 

ratio. 
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Decimal Order Input - X 1 - 18 

| 00 30K | 
Input a Sequence of Decimal Fractions - N 2 - 88 

00 10K 

1*0 192F 

50 L 

2k 30F 

kl F 

81 k¥ 

LO 11L 
3 32 6L 

iA- 11L 
k 50 F 

7k 11L 

S5 F 

kO F 

26 2L 

L5 F 

40 3F 

L5 7L 

iA 5L 

k6 7L 

LiX Li 

^0 L 

10 3^ 999F 
22 1L 

11 00 F 
00 10F 
26 ION 
00 10K 
00 F 
00 192F 
80 F 
00 F 
L5 3F 
kO 65F 



Store fractions 



Bring in parameters and convert 



LOCATION ORDER 

3 L5 ^F 

kO 66F 
k J 50 k~F 

7k kj 

5 S5 F 
10 IF 

6 iA L 
Jj-Q TF 

7 IfO 67F 
42 10L 

8 L5 ^F 
Lk k¥ 

9 itO 5F 
kl F 

10 L5 1L 
kO ( )F 

11 F5 F 
40 F 

12 LO k¥ 
36 999F 

13 L5 10L 
F^ 4f 

14 kO 10L 
00 30F 

15 26 10L 
26 2L 

26 IN 

00 13K 
00 63F 
26 190F 

00 52K 

I 
Square Root Routine - 21 



NOTES 
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Form 192 + n(n+l) / 2 



Generate unit matrix 



Read around delay 



LOCATION 



ORDER 



NOTES 
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1 00 22K 

22 00 IF 
00 IF 

23 00 Sk 
00 S^ 

2k 00 192F 

00 192F 

25 80 S5 
00 S5 

26 20 F 
00 F 

27 00 F 00 0100 
0000 0010 J 

28 80 F 
00 F 

29 00 F 
00 F 

30 JO S7 
lk S7 

00 76K 
Eigenvalue Roiifcine <*'~ h2 

1 1 

Change word 102 to read 

36 999F 

kO IF 
Change word 31 to read 

36 63F 

L5 6f 
Change word W- to read 

26 69F 

50 12QP 

00 69K 



LOCATION ORDER 

69 lb 10F 
LO 75F 

70 32 TIF 
LJ 10F 

71 22 120F 
L5 75F 

72 ^0 10F 
L5 T^F 

73 40 2F 

26 122F 
714- 7L i >-095F 

LL if095F 

75 3L ^095F 

LL V095F 

00 63K 
63 L3 6f 

32 61J-F 
6¥ 26 111F 

27 113F 

26 761 

00 82K 

Integer Print -*- Routine P k - 55 
Change word 22 to read 

82 4F 

92 6^3F 



NOTES 
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Tnese orders are necessary to avoid a 
division hang tip 



Deciinal point will now he printed out 



00 110K 

Fraction Print — Routine P 11 - 148 

Interlude II 
00 10K 
L5 65F 
^0 3F 
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10 



11 



12 



13 



1* 



15 



16 



IT 



L5 17L 
iA 66f 
4o to* 
L5 67F 
ij-0 5F 
F^ 20L 
kO 20L 
L7 (19^ )F 
L2 (192 )F 
32 7L 
L5 ^L 
00 20F 
k6 5L 
L5 20L 
F^ \L 
k-2 kL 
L5 20L 
LO k¥ 
32 3L 
19 38F 

Uo 6f 
F5 6f 
40 6f 
L5 5L 
^6 15L 
50 18L 
V 19L 
40 18L 
L2 ( )F 
32 11L 
26 999F 

/*yn -rn 

80 F 
00 F 



Bring out parameters 



Find smallest eigenvalue 



M 2 



By 13 



Find power of ten smaller than smallest 
eigenvalue 



■ — 
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18 


00 F 00 1000 
0000 0000 J 








19 


00 F 00 1000 








■*■ 


0000 0000 J 








20 


00 F 
00 IF 
26 ION 

00 10K 






i 

1 





kl F 
50 F 








1 


L5 (192 )F 


By k 


V 


j 




kO (192 )F 




Restore diagonal 
c ons ectit Ively 


elements X. 
1 


2 


F5 F 
^0 F 








3 


Kk F 
00 20F 








k 


Yk 1L 
kO 1L 








5 


L5 F 
LO ifF 








6 


36 1L 
L5 69L 








7 


^0 F 










L7 (193 )F 


By 11 s from 13 


Find smallest X. 




8 


L2 (192 )F 
32 10L 


By 10 






9 


L5 7L 
00 20F 








10 


hG 8L 
F5 7L 


From 8 








ij-2 7L 
F5 F 








12 

1 ' 


^4-0 F 

LO ^F 
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________ 

LOCATION , 


ORDER 
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15 


52 7L 
41 F 






14 


L5 8L 
46 18L 






15 


L7 70L 
10 IF 


From 18 




16 


40 70L 
F5 4% 




1 


17 


42 45L 




1 




L7 70L 




Find 2 smaller than smallest X. j 


18 


L2 ( )F 
56 15L 


By 14 


) 
I 

i' 


19 


L7 70L 


From 24 






66 (192 )F 


By 21 


2^ n / X. = d. 

1 1 1 j 


20 


S5 150F 




i 
1 




40 (192 )F 


By 22 


1 


21 


F5 19L 
42 1% 






22 


42 20L 
F5 F 






25 


40 F 
LO 4f 






2k 


56 19L 
4l 6f 






25 


41 9F 


From 67 






L5 6F 




Print row ntamber 


26 


JO 2F 

50 26L 






27 


26 82F 
92 151F 






28 


50 (192 )F 


By 50, from 54 






7J (S5) 


By 50 




29 


40 (150 )F 

19 18F 


By 52, fcy 65 


Evaluate d. a,. (1 = 1^ . .., n-l) 


50 


f4 28L 
4o 28L 







M 2 



LOCATION 



31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
4l 
42 

43 
44 

45 
46 

47 
48 



ORDER 



NOTES 
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19 18F 
L4 29L 
40 2% 
F5 9F 
40 9F 
LO 4F 
36 28L 
L5 6F 

40 7F 

41 F 
4l IF 
41 8F 

50 (130 )F 

7J (S5) 
50 IF 
74 69L 
L4 F 
40 F 
S5 192F 
40 IF 
19 18F 
F4 37L 
40 37L 
F5 8f 
4o 8f 
lo 4f 

36 37L 
L5 F 
50 IF 
00 ( )F 
40 10F 
S5 F 
40 11F 
L5 10F 
36 53L 
92 706F 



From 60 



By 42, "by 58, from 44 
By 42, Tsy 37 



By 17 



Perform matrix multiplication 



Count 



Scale up 



Print negative sign 



LOCATION 


ORDER | NOTES PAGE 9 


±9 


L3 11F 
36 53L 






50 


F5 10F 
1*0 10F 






51 


L5 11F 
L4 68L 






52 


4o 11F 
22 53L 






53 


92 963F 


From 48 


1 




L5 10F 


From 49 


Print integer j 


5^ 


JO s6 
50 5^L 




1 


55 


26 82F 








L5 11F 




Print fraction 


56 


JO S3 
50 56L 






57 


26 110F 
92 131F 






58 


L5 20L 
46 37L 






59 


F5 W 
40 7F 






60 


LO 4F 








32 35L 




Count 


61 


92 139F 
L5 40L 






62 


46 28L 
L5 20L 






63 


46 29L 
L5 5F 






64 


l4 4f 

40 5F 






65 


42 37L 
F5 6f 






66 


ko 6f 

LO kF 




1 
I 

i 



M 2 



LOCATION 


1— 

ORDER 




iroris 


PA.GE 10 


67 


36 25L 
26 l60F 




Count 




68 


80 F 
00 F 








69 


00 F 
00 IF 








70 


80 F 
00 F 

00 160K 











L5 79F 
kO F 








1 


L7 (192 )F 
L2 (195 )F 


By 5, from 7 
By 3 






2 


36 4L 
F5 1L 




Test for smallest d. 

1 


(d. =W\) 


3 


k2 1L 

26 6l 








k 


L5 IX 
00 20F 








5 


46 1L 
22 2L 








6 


F5 F 
40 F 


From 3 






7 


LO k¥ 
56 1L 




Count 




8 


L5 18F 
10 20F 




Bring out address of 


largest d. 


9 


k-2 l4L 
42 11L 








10 


L5 1L 
46 12L 








11 


46 l4L 
L7 ( )F 


By 9 






12 


L2 ( )F 
40 F 


3y 10 







M 2 



LOCATION 
13 

Ik 

15 
16 

IT 
18 



ORDER 



L3 F 
36 1TL 

L5 ( )F 
66 ( )F 
S5 11F 
50 15L 
26 110F 
22 17L 
92 66F 
92 7T0F 
92 979F 
OF F 

2b 10H 



NOTES 
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By 11 
By 9 



From 13 
From 16 



Test if \ g = ^ 



I 



Print out ratio 



Print out N 
Spaces 



